A mathematical model of turbulent drag reduction by high-molecular-weight 

polymeric additives in a shear flow 
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Drag reduction, or, what is the same, mean velocity increase in a turbulent flow at a fixed pressure 
drop through the addition of tiny amounts (several parts per million) of high molecular weight poly- 
mers (Thoms effect), is known already for more than sixty years. Rather long ago it was understood 
that this effect is related to supramolecular structures formed in the flow. Recent experiments 
by S. Chu, E.S.G. Shaqfeh and their associates, where the motion of supramolecular structures 
was directly observed, made it possible to understand and quantify the dynamic interaction of the 
polymeric structures with the solvent (water) flow. These results lead to the construction of a math- 
ematical model of the Thoms effect, based on the Kolmogorov(1942)-Prandtl(1945) semi-empirical 
theory of shear flow turbulence. This is the subject of the present Letter. 

PACS numbers: 47.10.A-, 47.85.Ib, 47.27.Ak, 47.57.Ng 
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Let us first introduce the basic equations of the present 
work. Turbulent shear flow is considered, so the mean 
velocity field is assumed to be 



Ui 



U{X2), U2 = U3 = 0. 



(1) 



The density of the solution is indistinguishable from 
the water density p, therefore under zero mean mass force 
the averaged momentum balance equation has the same 
form as in the absence of the mass force: 



Tl2 = 



(2) 



In these relations the coordinate xi is reckoned along 
the mean flow, the axis X2 is directed perpendicularly to 
the wall X2 — 0. Furthermore, T12 is the only non-zero 
component of the Reynolds stress tensor, r is the shear 
stress assumed to be constant across the flow, averaged 
viscous stress is neglected. The bars denote probabil- 
ity average (mean) values, the turbulent fluctuations are 
denoted by primes. The quantity it, — [j j p)^!"^ is the 
dynamic, or friction, velocity. 

The equation of the turbulent energy balance in the 
shear flow takes the form: 



dX2 



(3) 



Here, unlike in the customary derivations, the mass 
force is taken into account: F = 0, F' is the mass force 
fluctuation. Furthermore, e is the mean dissipation rate 
of the turbulent energy into heat: 



e = 2i/D'D' 



-{dau'p + dpu'^){dau'p + di3u'^,), (4) 



D is the symmetric part of the strain-rate tensor, v is the 
fluid kinematic viscosity, the summation over repeated 



Greek indices from 1 to 3 is assumed. The derivation of 
the equations (|ll2p follows the general lines of [ij. 

The turbulent flow is a cascade of vortices of various 
scales. The basic Kolmogorov's hypothesis [i] (cf. also 
the Prandtl paper Q) can be presented in the follow- 
ing form: at large Reynolds numbers the vortex cascade 
is self-similar. Hence, all dimensionless flow field prop- 
erties are universal. According to this hypothesis, all 
kinematic flow properties are determined by two kine- 
matic properties of different dimensions. Kolmogorov [3| 
and his direct followers (see [ij) selected as such quan- 
tities the local mean turbulent energy per unit mass 
6 = {u'l + U2 + u'^)l2 — and the local mean 

length scale of vortices, or the local external length scale 
of turbulent flow I, which is proportional to it. This 
hypothesis forms the basis of the (6, 1) semi-empirical 
model. The other possibility is (6, e) model which also 
came into wide use. (I used here the original Kol- 
mogorov's notation, often the turbulent energy is denoted 
by k). Furthermore, the momentum exchange coefficient 
K is introduced: 



K = Ti2/idu/dx2). 



(5) 



The relation ([5]) is not a new assumption. 
Dimensional analysis suggests the expressions for K 
and e (f3|, see also [H) via I and b: 



K = lVb, e = "f^^/^/l. (6) 

where 7 is a universal constant. The estimates in [l[ 
show, that 7 ~ 0.5. 

Thus we come to the following system of equations for 
the turbulent shear flow when mean mass force is equal 
to zero, but the mass force fluctuations are different from 



iVb- 



du 

dX2 



(7) 



- the mean momentum balance equation, and 

zV5(-^)2~74&!^+FTf^=0 (8) 

dX2 I 

- the mean turbulent energy balance equation. 

Once the basic equations have been presented we are 
in place to propose the model of Thoms phenomenon. 
Water with polymeric additives which has the property 
of turbulent drag reduction is not a genuine solution: 
supramolecular polymeric structures attaching the sol- 
vent molecules are formed in the mixture. The role of 
such structures in the Thoms phenomenon was suggested 
rather long ago [1, [1] . In the paper the photographs 
of wormlike structures in the carboxymethylcellulose (of 
rather low molecular weight ~ 70, 000) water solution 
were presented. An indirect confirmation of visco-elastic 
behavior of^polymeric supramolecular structures was ob- 
tained in [y]: a strong water jet directed on a metallic 
plate did not affect it for hours. Tiny polymer addi- 
tion leads to a fast abrasive-like wearing of the plate. 
However, only in the works by S. Chu, E.S.G. Shaqfeh 
and their associates the results of direct observations and 
quantitative experimental measurements of the motions 
and deformations of the supramolecular structures were 
obtained (see 0, S, 0] and the references presented in 
these papers). Chu, Shaqfeh and their colleagues ob- 
served complicated motions (" tumbling" ) of supramolec- 
ular structures, accompained by their deformation, over- 
turning, rotation, and adjustment of the elongated struc- 
tures to the local flow. They measured the characteristic 
time of tumbling, the "disentanglement time" 9. 

These experimental studies allow to suggest the basic 
hypothesis of the proposed model: 

The force F acting on the solvent (water) by the 
supramolecular polymeric structures, which consist of 
networks of polymeric molecules with attached solvent 
molecules, is proportional to the concentrations of the 
polymers. It is determined also by the instantaneous ve- 
locity fluctuation u' = u — u, and the disentanglement 
time 9. 

Dimensional analysis and symmetry considerations 
give 




FIG. 1; (color online) Dependence of X on different concen- 
trations C (see text for details). 

the density of solutions can be assumed to be equal to 
the density of the solvent, the first basic equation ^ 
remains the same as for pure solvent flow. 

However, the second basic equation is different. It 
takes the form: 

,^(|L,>-,«5I-?^.o (10) 

because, obviously, u'^ — u'^u'^ — 2b. 

This means physically that to the first two terms of 
the turbulent energy balance equation, representing cor- 
respondingly the inflow of the turbulent energy from 
the mean motion, and dissipation of the turbulent en- 
ergy into heat, a third term is added, which represents 
the work spent by turbulence for "disentanglement" of 
supramolecular structures. Now transform the equation 
(fTU| using dZl): 

Denoting by X , we obtain for X an algebraic 

equation of the fourth degree: 



-As— 



(9) 



the minus sign is due to the obvious fact that the 
force is a reaction to the velocity fluctuation. Also, A 
is a dimensionless constant, which could be, in princi- 
ple, included to the time 9. It was not done here, bear- 
ing in mind that the disentanglement time was indepen- 
dently introduced and measured by the experimentalists 
0,il. 

We summarize: for dilute solutions, when the interac- 
tion of supramolecular structures can be neglected, and 



X^ + 2CX^ -1 = 
where the dimensionless parameter C is: 

C= — -—. 

7^ u^,9 



(12) 



(13) 



If the polymer concentrations are zero, C — 0, and we 
return to the known result: X = 1; b = ul/^"^ = Const. 

For positive C, i.e. for a concentration different from 
zero, X < 1, b < u^/7^. This demonstrates qualitatively. 



without any assumption concerning Z, the reduction of 
the turbulent energy under the action of polymeric addi- 
tives. 

We come now to a quantitative model. The length 
scale I is an independent property, therefore the momen- 
tum exchange coefficient K — l^/b is also reduced un- 
der the action of polymeric additives. According to ([7]), 
du/dx2 = u1/K, i.e. the velocity gradient under con- 
stant shear stress (constant u*) is increasing. This leads 
to an increase of the mean velocity. 

Introduce into our consideration a specific expression 
for the length scale I. We accept, as it is done usually, 
I — Const ■ X2- According to classic approach Const — 
K7, where n is the von Karman constant, k — 0.42. More 
refined analysis (see [13, [HI [I2| ) gives for l{x2) a power 
law. At large Reynolds numbers, in the intermediate 
interval of X2 the function l{x2) can be represented as a 
linear function of X2 so that Const = Koo7, where Kqo = 
0.28. We emphasize that this distinction does not lead 
for the problem under consideration to any qualitative 
changes in analysis. So, we assume I — Koo^X2 so that 



C 



KqqAs X2 
7^ u^,9 



(14) 



The plot of X = 76^/^/M* as the function of C is pre- 
sented in Figure 1; it illustrates the decrease of the tur- 
bulent energy. For the momentum exchange coefficient 
K = iVb we obtain K — KooX2U*X{C), and the mean 
velocity is obtained by integration: 



dxo 



X2X{C) 



dC 



cx{cy 



(15) 



because C is, according to , proportional to X2 ■ It is 
possible to reduce the last integral to elementary func- 
tions. This leads to awkward relations, therefore it is 
worthwhile to consider the asymptotics, and to compute 
numerically the result in the intermediate range of C. In 
the case C << 1 (small polymer concentration, small 
a;2, or large disentanglement time 9) equation (|12p gives 
X{C) = l- C/2, so, according to ^ 



X2 



As 



u = In — - + X2 + Const. 



(16) 



In the case C >> 1 the first term in (fT^ is small 
in comparison with the second one, and we get X ~ 
(2C)-i/3. This leads to a power law velocity profile: 

u - Z{2K^As/-^^Y'^uJ^9-^'^xy^ + Const (17) 

For the intermediate values of the parameter C the sec- 
ond integral in can be represented as In C + G(C) + 
Const, where 




FIG. 2: Integral G(C) as a function of the concentration C. 



G{C) 



[l-^(g)] 
CX{C) 



dC 



(18) 



The plot of the function G(C) is presented at the Fig- 
ure 2. The formulae ([T5)) . allow us to calculate the 
velocity distribution at arbitrary values of C : 



■In 



\U*9y Koo ^ 



I^oqAs X2 

7^ u^9 



Const, (19) 



where the constant is determined by the boundary con- 
dition at the upper boundary of the viscous sublayer. The 
Const in equations (|16ll7ll9p depends on the Reynolds 
number {Re), and at Re ^ 00 it tends to —00. 

Now we will present a discussion and the conclusions 
of the current work. We constructed a semi-emprirical 
mathematical model of turbulent drag reduction by poly- 
meric additives. It is plausible that the same model 
will be adequate for shear flows of suspensions of small 
threads, where the phenomenon of turbulent drag reduc- 
tion was also observed. 

In a discussion of the preliminary version of the present 
work Professor P. Moin suggested that the decrease of 
the mean square vorticity w'^, which was also observed 
in the experimental studies can be explained along the 
same lines. This statement is correct. Indeed, the equa- 
tion for w'^ can be obtained in the same way as the 
equation for the turbulent energy balance: the opera- 
tion curl should be applied to the Navier-Stokes equa- 
tion, the obtained relation for the vorticity multiplied by 
2w = curlu, and averaged. The next step: the equation 
for the mean vorticity obtained by the operation curl ap- 
plied to the Reynolds equation (averaged Navier-Stokes 
equations) is multiplied by the mean vorticity, and the 
resulting equation subtracted from the previous one. In 
this argument instead of the mass force F the quantity 
I = curlF will enter. Obviously, similarly to the mean 



4 



value of the mean force F, the mean value I vanishes, 
due to symmetry. However, there appears an additional 
term in the equation for w'^ balance: w'l'. Due to the 
same reasons as in the case of mass force, the quantity I' 
should be proportional to sw' /d: 



where _B is a constant, analogous to A. Therefore in the 
equation for w'^ an additional term 

- Bsv?^/e, (21) 

will appear, analogous to the last term in equation pT|) . 
This demonstrates the reduction of w'^ under the action 
of the polymeric additives. 

The proposed mathematical model of turbulent drag 
reduction by high molecular weight polymeric additives 
was constructed under some simplifying assumptions, in 
particular, constancy of the concentration s, and con- 
stancy of the disentanglement time 6. Bearing in mind 
the degree of the accuracy of the Kolmogorov-Prandtl 
semi-empirical theory of the shear flow turbulence, these 
assumptions do not seem to be restrictive. 
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